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Abstract 

Superbursts are rare day-long Type I X-ray bursts due to carbon flashes on accreting neutron stars 
in low-mass X-ray binaries. They heat the neutron star envelope such that the burning of accreted 
hydrogen and helium becomes stable, and the common shorter X-ray bursts are quenched. Short 
bursts reappear only after the envelope cools down. We study multi-zone one-dimensional models 
of the neutron star envelope, in which we follow carbon burning during the superburst, and we 
include hydrogen and helium burning in the atmosphere above. We investigate both the case of a 
solar composition and a helium-rich atmosphere. This allows us to study for the first time a wide 
variety of thermonuclear burning behavior as well as the transitions between the different regimes in 
a self-consistent manner. For solar composition, burst quenching ends much sooner than previously 
expected. This is because of the complex interplay between the 3a, hot CNO, and CNO breakout 
reactions. Stable burning of hydrogen and helium transitions via marginally stable burning (mHz 
quasi-periodic oscillations) to less energetic bursts with short recurrence times. We find a short-lived 
bursting mode where weaker and stronger bursts alternate. Eventually the bursting behavior changes 
back to that of the pre-superburst bursts. Because of the scarcity of observations, this transition 
has not been directly detected after a superburst. Using the MINBAR burst catalog we identify the 
shortest upper limit on the quenching time for 4U f636-536, and derive further constraints on the 
time scale on which bursts return. 

Subject headings: accretion, accretion disks — methods: numerical — nuclear reactions, nucleosyn- 
thesis, abundances — stars: neutron — X-rays: binaries — X-rays: bursts 



1. INTRODUCTION 

Superbursts are day-long flares observed from neu- 
tron stars in low-mass X-ray b inaries (LMXBs; 
iCornelisse et al . 2000; Str ohmaver fc B rown 2002), that 
are attributed to the unstable thermonuclear bur ning 
of a carbon-rich la yer (|Cumming fc BildstenI 120011 : see 
iCooper et al.l [20091 for a discussion on alternative types 
of fuel). The thermonuclear runaway starts at a column 
depth oi y = lO^^gc m"^ to j7 = lO^ ^gcm"^, close to 



the outer crust (|Cumming et al.ll2006D . Because of the 
high temperature depend ence of carbon burning (e.g., 
iCaughlan fc Fowlej I1988D , the ignition is sensitive to 
crustal heating, and hence superburst observations place 
constraints on the combined crustal heating and neutrino 
cooling of the crust and core. For example, the super- 
burst from the classical transient X-ray source 4U 1608- 
522 requires the crust to h ave heated up fas ter than pre- 
dicted by current models (jKeek et al.ll2008l) . Because of 
the long typical recurrence time of approximately a year, 
superbursts are rare: so far 22 (candidates) have been 
det ected from 13 sources (e.g., [K eek & in 't Zand 2008: 
see'Kuulkcrs"2009'; Che nevez et a l. 2011; i n 't Zand et alj 
[2011: Serino et al..,2012clAsada et al...20Hl for recent dis- 
coveries; see also Table [2]). 

Ikeek@nscl.msu.edul 



The carbon fuel is expected to be produced by ther- 
monuclear burning of accreted hy drogen and /or helium 
highe r up in the atmosphere (e.g.. [Strohmaver fc Brownl 
I2002D . All known super bursting sources exhibit unsta- 
ble hydrogen/helium burning observed as short (up to 
^ 100 s duration) Type I X-ray bursts, as well as sta- 
ble burning, evidenced by high values of the a param- 
eter, the ratio of the integrated persist ent emission be- 
tween two bursts and the burst fluence (jin 't Zand et al.l 
|2003[ ). This burning takes place at a column depth of 
y ~ 10*gcm~^. Whereas stable hydrogen/helium burn- 
ing through the hot CNO cycle and the 3Q!-process pro- 
duce carbon, it is destroyed in bursts by CNO-cycle 
breakout reactions, followed by a series of helium cap- 
tures catalyzed by protons — the ap-process — and pro- 
ton captures with subsequent decays — the rp-proces s 
— creating heavier elements (jSchatz et al.l 1200 IL 120031 ) . 
Current models produce burst ashes with a c arbon mass 
fraction of typically 5% (jWooslev et al.|[2004D . whereas a 
slightly higher mass fraction is found for models that 
accrete material with a higher metallicity at a lower 
rate. Superburst observations, however, indicate a car- 
bon content that is closer t o 20% (Gumming et al. 2006). 
iMedin fc Gummingl ()2011f ) recently suggested that chem- 
ical separation by freeze out close to the outer crust may 
increase the carbon fraction. 
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If the superburst ignites in a sufficiently thick layer, 
carbon burning initi ally proceeds as a detonation 
(jWeinberg et al.l I2006D . This produces a shock that 
travels to the surface, pushing the envelope outward. 
The hydrogen/helium-rich atmosphere is heated by the 
shock (Weinberg fc Bildsten .2007,) and the subs equent 
fall-back of the outer layers (|Keek fc HegeHl20Il[) . This 
produces a bright precursor burst. Furthermore, the igni- 
tion conditions for a hydrogen/helium flash are reached. 
X-ray emission from this flash adds to the precursor 
burst, or it may account for the ent ire precursor emission , 
in case there is no strong shock (jKeek fc Hegeril2011f) . 
Precursor bursts have been identified in all cases when 
the start of the superburst w as observed an d the data 
were of sufficient quality (e.g.. iKuulker j[2"004f) . 

After the precursor burst, the envelope is sufficiently 
hot for all subsequently accreted hydrogen and helium 
to burn stably at the same ra te at which it is ac- 
creted: the bursts are quenched (ICornelisse et al.l 120001 : 
iCumming fc Bildsten|[200ll: iKuulkers et al.ll2002[ ). After 
the envelope has cooled down, unstable burning resumes, 
and Type I bursts return. Quenching times of one to sev- 
eral weeks are predicted (Gumming & Macbeth 2004). 
Because of non-contiguous observing schedules and the 
presence of frequent data gaps due to Earth occultations 
and South Atlantic Anomaly passages, the end of the 
quenching period and the reappearance of bursts have 
most likely not been observed. Only upper limits of more 
than one month have been derived for the quenching time 
(e.g. , iKuulkers 2004) . 

The return of hydrogen/helium bursts after a super- 
burst provides a unique opportunity to study the tran- 
sition of stable to unstable burning in the neutron star 
atmosphere. This transition is observed in a number of 
LMXBs, and is a ssociated with changes in the mass ac- 
cretion rate (e.g.. ICornelisse et al.ll20(33[ ). A higher rate 
implies a hotter envelope as well as a faster accumula- 
tion of accreted material, which leads to steady-state 
burning of hydrogen and helium. Current burst models 
predict that this transition occurs at a 10 times higher 
accretion rate tha n inferred from observations (e.g., 
iHeger et al.ll2007bh. alt hough increased crustal heating 
(e.g. JGupta et al.ll2007f) and rotati onally induced turbu- 
lent mixi ng |Piro fc Bildstenl l 2007D may reduce the dis- 
crepancy (jKeek et al.ll2009f r~AIternativelv. it may be the 
result of a higher local accretion rate. 

The transition of the thermonuclear burning behav- 
ior in the atmosphere after a superburst is solely due 
to the temperature change, whereas the mass accretion 
rate remains constant. Note that in practice the mass 
accretion rate can vary on the time scale of a superburst, 
which would provide an effect additional to the change 
in the burning behavior. Because the cooling after a su- 
perburst sets the thermal profile of the envelope, this 
poses strict constraints on the environment in which the 
first bursts ignite. Detecting the transition of stable to 
unstable burning after a superburst will allow us to de- 
termine the column depth and temperature at which the 
transition occurs. These are important ingredients for 
improving the current models. Unfortunately, the ob- 
servations of bursts after a superburst are scarce, and at 
the moment of writing the actual transition has not been 
observed. 

In this paper we provide one-dimensional multi-zone 



models of a neutron star envelope that undergoes a su- 
perburst. We create models with, respectively, a carbon- 
rich, a helium-rich , and a solar-composition atmosphere. 
We self-consistently simulate both the superburst and 
the burning processes in the atmosphere. The burst 
quenching is followed, as well as the transition to un- 
stable burning and the return of bursts. We focus on 
the nuclear reactions that are responsible for the differ- 
ent phenomena in the light curves. Finally, we combine 
a large set of observations to place constraints on burst 
quenching. 

2. NUMERICAL METHOD AND OBSERVATIONS 

The one-dimensional models of the neutron star enve- 
lope prese nted in thi s paper are similar to the superburst 
models by lKeek fc H eger (2011). We refer to that study 
for a detailed description of the code employed and the 
setup of the models. Here we describe the main prop- 
erties of our simulations , as well as the differ ences with 
respect to the models bv lKeek fc HegeJ (|201lD . Most no- 
table are the use of a different prescription for electron 
conductivity and a different implementation of accretion, 
which may cause small differences in the ignition condi- 
tions of superbursts. 

2.1. Stellar Evolution Code 

We create and evolve one-dimensional models of the 
neutron star envelope using the i mplicit hydrodynam - 
ics stellar evolution code KEPLER ([Weaver et al.l(l978f ). 
The version of KEPLER that we use is sim ilar to the 
version used in recent studies (e.g., Woosley e t al.l [20021 
12004 IHeger et al.ll2007a[) . We employ an adaptive one- 
dimensional Lagrangian grid in the radial direction, un- 
der the assumption of spherical symmetry. To follow 
the chemical evolution w e use a large adaptive nuclear 
networ k (R auscher et a ll 120031) including the hot-CNO 
cycle ([Wailace fc Wooslevf 1981h . the 3a-, rp-, and ap- 
processes, carbon fusion, and photodi sintegration . We 
take into account neutrino energy loss ([Itoh et al.| [l996). 
as well as radiative opacity and electron conductivity 
([Ibenlll975[ ). Mixing of the chemical composition by con- 
vection, semiconvection, and thermohaline circulation is 
implemente d as a diffusive process using mixing-length 
theory (e.g.. lClavtoT3ll968D . 

Mass accretion is implemented by increasing the pres- 
sure of t he outermost zone in the model at each time step 
(see also'Woosl ev fc Weavei<1984t[T"aam et all 1991 ). Pe- 
riodically a zone is added on top. Because of the details of 
the implementation, this causes a small dip in the light 
curve. We carefully check that this does not influence 
any of the conclusions that we draw about the bursting 
behavior. The light curves in this paper are corrected by 
omitting the brief time intervals when the dips occur. 

2.2. Model Setup 

We model the neutron star envelope above a radius of 
10 km and an enclosed mass of 1.4 Mq, up to the surface. 
The bottom of the surface zone is at a column depth of 
y = lO^gcm^^. This is well below the typical column 
depth for hydrogen/helium bursts of y ~ lO^gcm"^. 
The inner part of the model is formed by a 2 x 10^^ g 
iron substrate. It serves as a buffer into which heat gen- 
erated by the superburst can flow, and diffuse outward 
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on longer time scales. This ensures a correct late-time 
light curve. 

Crustal heating is taken into account by a constant 
luminosity at the inner boundary equivalent to Qb = 
0.2 MeVnucleon"^. Because of neutrino cooling in the 
substrate, the effective amount of crustal heating at 
the top of the substrate, close to the superburst igni- 
tion depth, is Qb = 0.1 MeVnucleon"^. This is in the 
range of typically assumed values for superbursters at 
accretion rates above 10% of the Eddington limit (e.g., 
Eaensel & Zdunik 1990; Gumming ct al. 2006). 

On top of the substrate we accrete a mixture of 80 % 
^^Fe and 20% ^^C. The latter is the typic al mass frac- 
tion of carbon that iCumming et al.l ()2006f) found from 
fits to observed light curves of hydrogen accreting su- 
perbursters, and ^^Fe is the most abun dant isotope in 
the a shes of hydrogen-rich bursts (e.g., iWooslev et al.l 
l2004f l. The ashes of helium-rich bursts consist mostly 
of somewhat lighter isotopes: ^^Si, ■^^S, '^^Ar, ■*°Ca (e.g., 
IJoss fc Lil[T980h . For our models, however, we use Fe, 
in order to have a consistent opacity of all envelope mod- 
els, giving comparable ignition conditions for the super- 
burst. 

Accretion of a hydrogen or helium-rich atmosphere is 
done only briefly before the superburst, because the com- 
putational expense of simulating the many short bursts 
that occur before the superburst is prohibitive. From a 
model where only a carbon-rich mixture was accreted we 
know already the moment of superburst ignition. 12.7 hr 
before this time, we replace the accretion composition 
by a solar (71% by mass ^H, 27% ''He, 2% "n), or a 
helium-rich mixture (98 % ''He, 2 % "N). This allows for 
several hydrogen and/or helium bursts to take place be- 
fore the superburst, ensuring equilibrium is reached for 
the effects of chemical inertia ([Wooslev et al.|[2004i) . The 
change of accretion composition does not affect the igni- 
tion of the superburst, because it depends on the pres- 
sure at the bottom of the carbon layer, which continues 
to increase at the some rate. Moreover, at that time the 
initial phase of the carbon runaway has already started. 

An accretion rate of M — 5.25 x 10~^ yr~^ is used. 
For an atmosphere of solar composition on a neutron star 
of 1.4 M© this corresponds to 30% of the Eddington- 
limited rate AfEdd = 1-75 x 10"* M© yr~^ 

The presented results are not corrected for the redshift 
due t o the neutron star's gravity (see also lKeek fc Hegeil 
1201 It ). Our Newtonian model has the same surface grav- 
ity as when general relativity is taken into account for a 
star with the same mass and a 11.2 km ra dius, which has 
a gra vitational redshift of z — 0.26 (e.g.. IWooslev et al.l 
[200l . 

Based on the atmosphere composition, we refer to the 
simulations as model "C", "He", and "H" (Table [1]). 

2.3. MINBAR Catalog of Observations 

To derive observational constraints on the bursting be- 
havior of superbursters we employ version 0.51 of the 
Mult i-INstrument Burst ARchive (MINB AR; |Keek et al.l 
120101 : see http: / / users.monash.edu.au/ ~dgallow/minbar| 
for more details) . i'his catalog contains the results ot the 
analysis of 4,192 observed Type I X-ray bursts from 72 
sources as well as 27,340 pointings on 84 sources. The 
observations have been performed with the Wide-Field 



Cameras (WFCs) on bo ard the Beppo Satellite per As- 
tronomia X ( j?eppo5'^X: iCornelisse et al.ll2003f ) and the 
Proportional Counter Array (PC A) on board the Rossi 
X-Ray Timing Explorer (RXTE: ICallowav et al.l [2008l. 
Both instruments are sensitive in a similar energy range 
above 2keV. Because of its larger collectin g area, the 
PCA is more sensitive to faint bursts (e.g., iKeek et "ahl 
l2010h . 

MINBAR comprises the largest collection of X-ray 
bursts available, but we refer to the literature if more 
constraining observations have been reported. 

3. RESULTS 

3.1. Hydrogen/Helium Atmosphere Models 

Initially the accretion of carbon-rich material is simu- 
lated to self-consistently build up a neutron star envelope 
close to the ignition of a superburst. Approximately half 
a day before the thermonuclear runaway of carbon burn- 
ing sets in, the accretion composition is changed. We 
create three different models: one accreting solar compo- 
sition (model H), one with helium-rich material (model 
He), and for comparison one where we retain the carbon- 
rich atmosphere (model C). From that moment we fol- 
low both the carbon burning in the ocean and the hy- 
drogen/helium burning in the atmosphere. Several nor- 
mal X-ray bursts occur before the superburst (Fig. [TJ 
[2]). During the superburst decay, hydrogen and helium 
burn stably: bursts are quenched. Once the envelope 
has cooled sufficiently, unstable burning resumes. At the 
transition, oscillatory burning takes place (marginally 
stable burning). The first bursts after the superburst 
are less energetic and have shorter recurrence times than 
the bursts before the superburst. When the envelope 
cools down further, the recurrence times become longer 
and the bursts become as energetic as they were before 
the superburst. For the models H and He we calculate 
over 900 bursts per model. 

3.2. Superburst 

The superburst occurs 1.28 years after the start of 
the simulations, after accreting a column of y = 1.1 x 
lO'^^gcm^^. Its ignition is not in-phase with the occur- 
rence of hydrogen/helium bursts (Fig. [U [5]): flashes in 
the envelope do not trigger the runaway carbon burn- 
ing. The superburst occurs slightly earlier in the models 
H and He compared to C: 1.1 minutes for helium accre- 
tion and 5.2 minutes for solar composition. This may 
be due to heating from hydrogen and helium burning, 
but the differences in recurrence time, and therefore igni- 
tion depth, are too small to cause substantial differences 
in, for example, the superburst fluence or other observ- 
ables. The amount of time before the superburst during 
which we simulate the hydrogen/helium bursts may not 
have been long enough to bring the model into thermal 
equilibrium, and the deeper layers could still be increas- 
ing somewhat in temperature because of heating by the 
bursts. Because most of the energy produced by the 
bursts is radiated away from the surface, however, this 
introduces only a minor deviation in the temperature at 
the carbon ignition depth, which can be modeled by a 
slightly higher effective Qb- 

The superbursts start with a detonation, which drives 
a shock toward the surface, and produces a brief shock 



4 



Keek, Heger, & In 't Zand 




-8 -7 -6 -5 -4 -3-2-10^ 5 10 15 20 24 25 26 27 28 29 30 31 32 

Time since start superburst (hours) 



Fig. 1. — Light curve of models H and C. Breaks on the axis indicate a change in time scale. Only a small fraction of the simulated 
bursts after the superburst are shown. 
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Fig. 2. — Light curve of models He and C. Breaks on the axis indicate a change in time scale. Only a small fraction of the simulated 
bursts after the superburst are shown. 
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Fig. 3. — Logarithmic light curve of superbursts from models He, H, and C. The curves end when the burst-quenching period ends, or, 
in the case of a carbon-rich atmosphere, when the next superburst occurs. 
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TABLE 1 

Properties of Superbursts and Subsequent Bursts for 
Different Atmosphere Compositions 



Model 


c 


He 


H 


Accretion composition 


0.2 ^^C 


0.98 -^He 


0.71 


(mass fractions) 


0.8 ^''Pe 


0.02 i^N 


0.27 "'He 








0.02 


ipcak(103^ ergs-i) 


8.7 


8.8 


9.3 


Eshockbr (10^^ erg) 


7.7 


4.4 


1.4 


i-shockbr (10^'' ergs-1) 


10 


6.2 


1.3 


Precursor duration (s) 


1.2 


3.2 


5.7 


^precursor (lO'^'^ ^^S) 


0.75 


1.5 


1.3 


^minimum (®) 


20 


25 


42 


*quench (days) 




11.3 


1.1 


Pose (minutes) ^ 




20 


5.0 


treturn (days) 




115 


35 



Time of minimum luminosity between precursor and peak. 
^ Period of oscillatory burning. 
Time when burst properties return to pre-superburst values. 



breakout peak in the hght curve (Fig. [3|), followed on a 
dynamical time-scale by a precursor burst. The super- 
burst peak is reached 6.4 minutes after the superburst on- 
set. Stable hydrogen (helium) burning produces a peak 
luminosity, Lpeak, that is 5 % (0.4%) higher for model 
H (He) compared to C (Table [l]). The decays follow the 
same two-component power-law profile, with again hy- 
drogen burning raising the luminosity by a few percent. 
The light curves for the models H and C deviate from the 
cooling curve of C when the luminosity becomes compa- 
rable to that from stable hydrogen or helium burning. 
The hydrogen or helium burning contribution to the to- 
tal luminosity exceeds 10 % after roughly ~ 2 x 10^ s for 
model H, and after ^ 1 x 10^ s for He. For all models 
a small amount of carbon burns at the bottom of the 
newly accreted material. The burning time scale, how- 
ever, is much longer than the accretion time, so this is 
not steady-state burning, but leads up to the next super- 
burst. 

3.3. Shock Breakout 

After the superburst detonation, a shock travels out- 
ward from the bottom of the carbon-rich layer. Once it 
reaches the surface, a shock breakout peak is produced 
in the light curve (Fig. [3]). We determine its maximum 
luminosity, Lghockbr, and its fluence, i?shockbr, as mea- 
sured within 5 x 10~^s after the onset (Table [1]). The 
differences in Lshockbr trace variations in the opacity of 
the outer atmosphere. 

The shock accelerates the outer zones, where the den- 
sity is lowest, and the shock over-pressure is highest, to 
a substantial fraction of the speed of light within a very 
short time interval. This behavior is likely not well re- 
solved by our model, and at times introduces large vari- 
ability in, for example, the luminosity in the outer zones 
during the shock breakout and the start of the subsequent 
precursor burst. During this part of the superburst we 
take the mean luminosity of the 10 outermost zones to 
reduce the effect on the light curve, although some of the 
introduced variability remains visible (Fig. [3]). 

3.4. Precursor 

Most of the energy of the shock is used to expand the 
outer layers. They fall back on a dynamical time scale of 
10~^ s, and dissipate that energy into heat, which pow- 
ers a precursor burst that reaches the Eddington limit. 




5.0x10"' 1.0x10"* 1.5x10"* 2.0x10"* 

Time since start superburst (seconds) 



Fig. 4. — For model He: energy generation/loss (color scale) in 
the neutron star envelope as a function of time in a short interval 
around the superburst onset. Green hatching indicates convection 
(Rayleigh— Taylor instabilities). 
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Fig. 5. — Same as Fig. |4]for model H. 

At the start of the precursor burst, the fallback causes 
the material to settle while undergoing a damped oscil- 
lation. This creates corresponding variability of super- 
Eddington luminosity as high as 3.7x 10'*''ergs~^ (Fig.|H 

The precursor light curves are dependent on the atmo- 
spheric composition (Fig. [3]). The heating of the atmo- 
sphere instigates the burning of hydrogen and/or helium 
upon fallback (Fig. 31 [SI), leading to a thermonuclear run- 
away. This adds to the precursor fiuence i?prccursor (Ta- 
ble [1]). In model H 26% of the 7.4 x 10^ g column of 
solar composition is burned, and in model He 67 % of 
the 9.0 X 10^ g helium column. Because of the difference 
in energy yield of the hydrogen and helium burning nu- 
clear reactions, the total energy released and, hence, the 
fluence is very close for these particular compositions. 

The luminosity reaches a plateau caused by photo- 
spheric radius expansion (PRE, Fig. [S]). We compute 
the duration of the PRE phase as the time from the pre- 
cursor onset to the time when the luminosity drops below 
90% of the plateau value (Table [T]). 

After the PRE phase, the luminosity drops quickly, 
reaching a minimum at tminimum (Table [1]), and sub- 
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sequently climbing to the superburst peak at approxi- 
mately 400 s. 

The precursor light curve contains some irregularities. 
In the figures we have filtered out strong variations oc- 
curring from one time step to the next, that are due to 
the outer few grid points, and not representative of the 
overall behavior of the model. 

3.5. Bursts in a Helium-rich Atmosphere 

First we discuss the effect of the superburst on the ther- 
monuclear processes in model He. After the superburst 
precursor, helium burning in the atmosphere proceeds in 
a stable manner. Bursts only reappear once the envelope 
has cooled down sufficiently after 11.3 days (Fig. [2]). 

Before bursts resume, the burning is marginally sta- 
ble when the temperature in the atmosphere drops to 
3 X 10^ K, leading to oscillations in the light curve with 
a period close to 20 minutes (Fig. The oscillations 
have a small amplitude of at most approximately 10 % 
of the average luminosity, and numerical noise damp- 
ens the oscillations during certain periods. Once the 
bursts start, the burst peak flux increases, and the third 
burst already has a peak luminosity comparable to the 
bursts before the superburst (Fig. [HI [7]). During 2.3 days 
we find both bright bursts with a peak luminosity of 
Lpoak — 4 X 10'^® erg s"""^ (the Eddington luminosity for a 
hydrogen-deficient atmosphere) and weaker bursts with 
Lpoak — 3 X lO^^ergs"^. The bright as well as the weak 
bursts have a relatively slow rise of the light curve. Dur- 
ing this time 3a is the predominant nuclear process, pro- 
ducing ^^C. In the bright bursts 3a raises the temper- 
ature sufficiently for a-captures to take over, causing a 
faster rise of the luminosity, and producing mostly ^^Si. 
The bright bursts heat the envelope sufficiently for ad- 
ditional stable helium burning to take place, followed by 
a weak burst. Stable helium burning reduces the he- 
lium content of the atmosphere, leading to a burst with 
a lower Lpcak- The weak burst does not heat the en- 
velope enough for additional stable burning, such that 
the next burst is again bright, and ignites after a shorter 
recurrence time. While the envelope continues to cool 
from the superburst, the number of bright bursts in be- 
tween weak bursts increases, until after 2.3 days the weak 
bursts disappear. 

Comparing the profiles of individual bursts after the 
quenching period, the weak bursts have a relatively long 
rise of ~ 5s (Fig. [8]). The early bright bursts share this 
slow rise to a similar luminosity as the peak of the weak 
bursts, but then transition in a fast sub-second rise to 
the Eddington luminosity. Later bright bursts show an 
initial "bump" that decreases in duration over time. The 
bursts from before the superburst show a slow rise com- 
ponent with a duration of 0.4 s. The decay of the bursts 
becomes longer as the recurrence time increases, refiect- 
ing the longer thermal time scales of increasing ignition 
depths. 

The flux from the accretion process at a rate of M = 
5.25 X 10~^MQyr~^ for a neutron star of 1.4M0 with 
a 10km radius is 4.35 x lO^^ergs"^ (assuming isotropic 
emission and a 100% efficient accretion process). The 
initial oscillations and weak bursts have a lower Lpcak, 
but the brighter bursts outshine the accretion flux by 
approximately a factor 10. 



The superburst burns most carbon out to y ~ 7 x 
lO^gcm"^. Between this depth and the ignition depth 
of the helium bursts at y ~ 1 x lO^gcm"^, a car- 
bon mass fraction of 0.1 survives the runaway burn- 
ing, but burns on a longer time scale. Carbon burns 
through ^^C(^^C, Q;)^'^Ne and subsequent a-capture re- 
actions, producing predominantly ^^Si. When we stop 
the simulation, the carbon mass fraction of this material 
varies from 2 x 10"'^ down to 2 x 10~^. 

Directly after the superburst, during the burst quench- 
ing period, helium burns stably to carbon by the 3a- 
process. Captures of a on ^^C and ^'^N from the accretion 
composition produce a limited fraction of ^^O, ^^F, and 
heavier isotopes up to magnesium (Fig.|9]). Over time, as 
the atmosphere cools down, the a-capture rates reduce, 
and only lighter elements up to neon are produced. The 
main product of nuclear burning during the quenching 
period, however, is ^^C, with a mass fraction of 95 %. Af- 
ter the quenching period, this material is compressed to 
higher densities, and ^^C(^^C, a)^"Ne reduces the carbon 
mass fraction. Subsequent a-capture reactions produce 
mainly ^^Si, '^^S, and '^^Ar. The next superburst ignition 
occurs close to the bottom of this layer, and depends on 
the remaining carbon fraction. When we end the simu- 
lation, the ^^C fraction is 23%, and still dropping. 

During the quenching period, most burning takes place 
close to y ~ lO^gcm^^, but not all helium is burned. 
Some survives down to a depth of y ~ 10^"gcm~^. The 
low mass fraction {Y < 10^^) makes 3a inefficient, but 
series of a-captures {^^C{a,j)^^0{a,j)^°Ne{a,j)^^Mg) 
still occur. After stable burning at ?/ ~ 10® gcm~^ ends, 
the a-captures in the deeper layers continue until he- 
lium is fully depleted. During this time, even though the 
helium mass fraction is low, because of the large total 
mass down to y ~ 10^"gcm~^, the captures contribute 
a substantial fraction of the total generated energy at 
any given time. This effectively slows down the cool- 
ing for up to 8.1 days, keeping the recurrence time of 
the bursts constant. Once helium is depleted at these 
depths, the cooling continues, and the burst recurrence 
time increases over time. 46.3 days after burst resump- 
tion the recurrence time is 88 minutes, whereas the pre- 
superburst recurrence time was 100 minutes. We stop the 
simulation here, but extrapolating the trend, the original 
recurrence time will be recovered approximately 115 days 
after burst resumption. 

The model produces bursts with a spread of 50 % in 
ipcak (Fig-El) with a quasi-periodic behavior. To a lesser 
extend these variations can also be seen in irccur- This 
is caused by the dependence of the ignition conditions 
of a burst on the previous bursts, based on, e.g., the 
fraction of helium that was burned, the heat deposition, 
and compositional inertia. 

3.6. Bursts in a Solar-composition Atmosphere 

Burst quenching is also present in model H (Fig. [1] 
I10[) . Directly following the superburst, hydrogen and he- 
lium burning is stable. Helium burns by the 3a-process, 
providing CNO elements which facilitate hydrogen burn- 
ing by the hot-CNO cycle. Break-out from this cycle 
by ^^0(a,7)^^Ne allows for a rapid series of proton cap- 
tures and /3-decays (rp-process) to produce mostly ^^Ga 
and ^^Ge, whereas a large number of isotopes with mass 
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Fig. 6. — Light curve of model He at the time when quenching stops and bursting resumes. The dotted Hne indicates the level of the 
accretion flux. 
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Fig. 7. — Peak luminosity Lpeak and recurrence time trecur of 
bursts after the quenching period for model He. The dotted lines 
indicate the mean values for bursts before the superburst and the 
luminosity of the accretion process. 

numbers in the range of 59-72 contribute mass fractions 
of several percents (Fig. [TT|) . 

Similarly to model He, after the superburst a layer 
with a ^^C mass fraction of 10% remains between y ~ 
1 X lO^gcm"^ and y ~ 7 x lO^gcm"^. It burns on a 
longer time scale down to a mass fraction of 2 x 10"'*. Hy- 
drogen and helium burning during the burst quenching 
period produce carbon at a mass fraction of 5 %. Nuclear 
burning on longer time scales reduces the carbon mass 
fraction in this layer to 10~^. 

Burst quenching is much shorter than in the model 
with a helium-rich atmosphere: bursts resume after 
1.1 days. As the superburst cools, the temperature 
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Fig. 8. — Comparison of light curves of helium bursts aligned 
on the peak. The smallest burst occurred 34 minutes after bursts 
resumed, and the one with the bump before the peak occurred 
after 2.3 days. The burst with the longest decay occurred before 
the superburst. Superburst emission has been subtracted for all 
bursts. 
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quenching period, down to the bottom of the layer of ashes from 
stable helium burning. Only the most abundant isotopes are 
shown. 



drops below 4.3 x 10^ K, and the breakout reaction 
^^0{a, 7)*^Ne becomes less efhcient than the 3a-process. 
Whereas at higher temperatures the break-out reactions 
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Fig. 10. — Light curve of model H at the time when quenching stops and bursting resumes. The dotted hne indicates the level of the 
accretion flux. 
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Fig. 11. — Composition of the envelope at the end of the burst- 
quenching period, down to the bottom of the layer of ashes from 
stable burning of solar-composition material. Only a selection of 
the most abundant isotopes is shown. 

quickly removed CNO elements, now the CNO mass 
fraction is growing. This causes an increase in the he- 
lium production through the CNO cycle, and results in 
the increase of the energy generation rate of both the 
CNO cycle and the 3a-process. The involved reactions 
raise the atmosphere temperature until once more the 
^^0(q!, 7)^^Ne break-out can efficiently remove most of 
the CNO elements, providing seed nuclei for the rp- 
process, which captures most of the hydrogen. With 
hydrogen and the CNO isotopes gone, the CNO cycle 
and rp-process switch off, reducing the helium produc- 
tion, and thereby the 3a rate, allowing the atmosphere 
to cool down. Note that while hydrogen is depleted, he- 
lium is not. As fresh hydrogen and helium are accreted 
and mixed in from layers closer to the surface, the cycle 
repeats itself. This produces a series of oscillations in the 
light curve, that announces the end of the steady-state 
burning of the burst quenching phase (Fig. fTU)) . The time 
scale for the oscillations is 5.0 min. 

The oscillations increase in amplitude over time. Af- 
ter 2.1 hr the ^^0(a, 7)^^Ne breakout and subsequent rp- 
process increase the temperature during the oscillations 
sufficiently for the ^'*0(a,p)^^F(p, 7)^*Ne breakout reac- 
tions initiate the ap-process, which causes a faster rise 
of the luminosity, producing small bursts instead of os- 
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Fig. 12. — Comparison of light curves of hydrogen/helium bursts 
aligned on the peak. From small peak to high, the bursts occurred 
1.5 hr, 4.6 hr, 2.7 days, and 13.1 days after bursting resumed, 
respectively. The brightest burst occurred before the superburst. 
Superburst luminosity has been subtracted for all bursts. 

dilations. For ^ 1.5 hr oscillations and small bursts al- 
ternate, until an equilibrium is reached and only small 
bursts occur (Fig. fTU]) . The bursts following an oscillation 
have a shorter recurrence time. Initially the bursts have 
a slow rise, similar to the oscillations, and a faster decay. 
Over time, the rise shortens, and the decay lengthens 

(Fig.ini). 

When the oscillations have disappeared, and only 
bursts remain, the recurrence time is constant at 
8.2 minutes until approximately 1.2 days after burst re- 
sumption. As in model He, this is due to the burning of 
residual helium at greater depths, which delays the cool- 
ing of the atmosphere. Although the recurrence time is 
constant, the peak luminosity Lpcak increases with time. 

Afterward both the recurrence time and ipeak increase 
with time (Fig. After 13 days trccur — 39 minutes 

is reached. The burst peak luminosity is then at 77% 
of the pre-superburst level. The burst light curve ap- 
proaches the profiles of the bursts before the superburst. 
At this time we end the simulation. Extrapolating, the 
recurrence time of 52 minutes of the bursts before the 
superburst will be reached approximately 35 days after 
burst resumption. 

For 1 X 10^ 5, i 5, 4 X 10^ s there are some 'oscilla- 
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Fig. 13. — Peak luminosity ^^pcak and recurrence time trccur of 
bursts after the quenching period for model H. The dotted lines 
indicate the mean values for bursts before the superburst and the 
luminosity of the accretion process. 

tory' variations in both tiocur and Lpcak, which are hkely 
of similar origin as the late-time variations in the burst 
properties in the helium-rich atmosphere (Fig. I13l) . 

3.7. Observational Limits on Burst Quenching 

The first normal burst observed from a source after a 
superburst provides an upper limit to the burst quench- 
ing time. For all known (candidate) superbursts, we 
identify the first detected burst either from MINBAR, 
or from the literature (Table [5]) . The superbursts from 
4U 1735-444 and 4U 1820-303 occurred when (or close 
to when) the source exhibited a persistent flux where no 
or very few bursts have been detected for the respective 
sources. Most probably hydrogen and helium were un- 
dergoing steady-state burning already before the super- 
burst started. For this reason we exclude these super- 
bursts when we investigate constraints on burst quench- 
ing. Furthermore, we do not consider the recent su- 
perbursts from EXO 1745-248, SAX J1747.0-2853, and 
SAX J1828. 5-1037, which are still being analyzed at the 
moment of writing, nor those from GX 17+2, which has 
an at ypically high mass accretion rate (|in 't Zand et all 
[200l . 

Two superbursters are so-called ultracomp act X-ray bi- 
naries (UCXBs; e.g.. lin 't Zand et alll2007D . where the 
accreted material is thought to be hyd rogen deficient, 
but m ay contain helium: 4 U 0614+0 91 (Kuulkers et al.' 
20101) and 4U 1820- 303 (iStrohmaver fc Brow n 2002; 



in 't Zand erani20lH see also iCummi^ |2003( 1 . From 



4U 0614+091 two bursts were observed in the period af- 
ter the superburst: af ter 18.6 and 34.9 days, respectively 
(|Kuulkers et al.l [20101 ). Both bursts reach a bolometric 
peak flux that is consistent within la. The first burst is 



TABLE 2 

Observational Limits on Burst Quenching for All Known 
(Candidate) Superbursts 



Source 


Time (MJD)'' 


First burst 
(d)b 


tcxp (hr)'^ 


4U 0614+091 


53441.70 


(-0.26) 


18.6'* 


0.9 


4U 1254-690 


51187.39 


124.7 


25.5 


4U 1608-522' 


53495.08 




99.8<= 


2.6 


4U 1636-536 


50253.61 


(-0.07) 


96.3 


0.0 




50642.37 


(-0.07) 


68.6 


2.1 




51324.21 


(-0.06) 


15.0 


33.9 




51962.70 


22.8 


9.7 


KS 1731-260 


50349.42 




34.4 


48.0 


4U 1735-444 


50318.13 


(-0.02) 


374.3* 


130.6 


EXO 1745-248* 


55858.53 


(-0.06) 


- 


- 


GX 3+1 


50973.04 


(-0.08) 


94.2 


2.5 


SAX J1747.0-2853' 


55605.54 




25.1* 




GX 17+2 


50340.30 


(-0.03) 


2.2 


36.6 




51444.10 


(-0.04) 


_g 


132.0 




51452.33 


2.3 


113.4 




51795.34 




12.8 


100.5 


4U 1820-303 


51430.07 




167.0 


113.7 




55272.72 


(-0.06) 






SAX J1828. 5-1037' 


55877.34 


(-0.06) 






Ser X-1 


50507.08 


(-0.03) 


34.6 


45.9 




51399.14 


(-0.07) 


309.4 


1.5 




54753.28 


(-0.07) 


54.7 


0.0 



Time of the observed start of the superburst (see Sect. [T] for 
references). In case this directly follows a data gap, the duration 
of the gap is indicated in parenthesis. 

^ Using bursts from MINBAR, unless indicated otherwise. 

Total exposure time in first 35 d after the superburst; "-" indi- 
cates the superburst occurred outside of the period covered by the 
catalog. 

'*IKuulkers et al.l ||20T0|). 

° "Keek et al. ( 2008). 
' Linares et ah, 12011a ). 

s There were no bursts observed before the subsequent superburst. 
* Superburst took place when hydrogen and helium were likely 
burning in steady state. 

' Transient source with outburst durations shorter than the super- 
burst recurrence time. 



more symmetric with a rise of 5 s and a decay of 2.0 s, 
whereas the second one has a faster rise of 1 s and a 
slower decay of 13.0 s. It should be noted, however, that 
the rise and decay time scales were derived from differ- 
ent energy bands, with the first burst being observed at 
higher energies, causing part of the shorter decay time. 

The other superbursting sources are thought to accrete 
hydrogen-rich material. As the model with the hydrogen- 
rich envelope predicts normal bursting behavior to re- 
turn in 35 days, we identify the bursts from MINBAR 
observed in that period, as well as the net exposure time, 
^cxp, of all observations, including those without bursts 
(Fig. [T31 Table [2). The earliest observed normal burst 
occurred 15.0 days after the superburst from 4U 1636- 
536 on MJD 51324.21. After most superbursts the total 
exposure time is short, and observations are infrequent. 
Therefore, we combine the data for these sources in 1 day 
time bins to be able to distinguish changes in the burst 
behavior on similar time scales as exhibited by our mod- 
els. To place constraints on the burst recurrence time, 
^rccur, we usc the combined tcxp and number of bursts 
during each day, taking into account the presence of data 
gaps by Earth occultations and by passages through the 
South Atlantic Anomaly. 

Data gaps reduce the fraction of the observation time 
iobs when bursts can be observed: icxp = V^ohs, with rj the 
observation efficiency. For rj < 1 bursts can be missed. 
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Fig. 14. — Cumulative net exposure time during 35 days after a superburst for hydrogen-accreting sources. We exclude sources with 
less than 3hr of exposure time, and those for which hydrogen/helium burning is stable, as well as GX 17+2 (see also Table [2l . For each 
superburst we indicate the source name and the time of the superburst (MJD). Circles indicate the detection of normal bursts. 

To determine the expected trocur before the super- 
bursts, we search MINBAR for the shortest time interval 
between two subsequent bursts that occurred at a similar 
level of persistent flux as the superburst. For six super- 
bursts we find time intervals between 1.4 hr and 3.6 hr, 
whereas for others we find time intervals of 8.0 hr and 
longer. The longer trecur are most likely because of a 
short total exposure time on a given source at the level 
of persistent flux of interest, and not representative of 
the actual trccur- Indeed, the shorter time intervals are 
found for most of the superbursts in Fig. [TH which have 
the longest cumulative exposure times. Therefore, we 
take 1.4-3.6 hr to be the range of expected values for 
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Fig. 15. — Observational limits on the burst recurrence time 
trccur from Combined observations after the superbursts from 
Fig. 1141 in 1 day time bins. The region between the gray areas 
is the 90% confidence region for trccur. The solid line is the net 
exposure time t^xp per day, and the hatched region indicates the 
range of minimum values of trecur observed for bursts at similar 
persistent flux as at the time of the superbursts. 

and when N bursts occur, n < N are detected. Using 
a binomial distribution, the probability P of detecting n 
bursts out of N is 



P{n;N,rj) 



N 



N-n 



Since n and 77 are known for each observation, we use P 
to identify those values of N for which the probability 
of detecting n is largest. After normalization, we deter- 
mine the 90% confidence region for N, which we use to 
constrain ^rccur = ^obs/A^- 

As additional constraint, if no bursts are observed, we 
require that irccur 

is longer than the uninterrupted part 
of a pointing. For a 96 minute satellite orbit we esti- 
mate this as (96 minutes) x ry. Furthermore, in our su- 
perburst selection only for the superburst on 51324.21 
from 4U 1636-536 more than one burst was observed in 
the following month. The closest pair was separated by 
17.8 hr. We use this as upper limit for irocur at 24 days 
after the superburst onset. The combined constraints are 
presented in Fig. [T5l 



^rccur- 

Comparing the expected pre-superburst trecur to the 
limits on t^ccur (Fig. [T5)) . we see that trocur < 5.5 hr is 
strongly disfavored during the first 4 days, with trccur 
likely exceeding 12.3 hr during the first day. Therefore, 
these superbursts caused a substantial increase in the 
burst recurrence time, or possibly quenched bursts alto- 
gether. Furthermore, for most 1 day time bins the lower 
limit for trocur is at least 0.5-1.0 hr, which disfavors bursts 
with shorter trccur- 

The strongest constraints are due to a burst from 
4U 1636-536 on day 26: 0.5 hr < Wur < 1-5 hr. This 
suggests trccur IS back at pre-superburst values after 26 
days. Only 5 days earlier, the lower limit was irecur > 
4.2 hr, and on day 34 ^rocur > 3.2 hr: both at the up- 
per part of the pre-superburst range. Therefore, the al- 
lowed values for irccur within the pre-superburst range 
vary somewhat with time, which may be due to differ- 
ences in the persistent flux from one superburst to the 
next as well as variations during the month after a given 
superburst. 

The majority of the bursts indicated in Fig. [T3] followed 
the superburst on MJD 51324.21 from 4U 1636-536. All 
those bursts have similar properties, and exhibit PRE, 
indicating that the Eddington limit was reached. The 
persistent flux was approximately 6 % of the burst peak 
flux. The burst properties are typical for bursts from this 
source at that level of persistent flux. Burst properties 
not changing between days 15 and 26 suggest that normal 
bursting behavior was resumed at most two weeks after 
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the superburst. 

4. DISCUSSION 

We present one-dimensional muhi-zone simulations of 
the neutron star envelope, where we study the effect of a 
superburst on a helium-rich and on a solar-composition 
atmosphere, where matter is accreted at a rate of M = 
5.25 X 10~^ M pyr~^. The simulatio ns are similar to those 
presented by iKeek &: Hegeii ()201in where a carbon-rich 
layer is build up by accretion until a superburst ignites. 
In the present work we replace the accretion composition 
shortly before the superburst to build an atmosphere of 
either a helium-rich or a solar composition. The results 
are summarized in Table [TJ 

We find that the hydrogen/helium-rich atmosphere 
changes the shock breakout and precursor peak height 
and duration, whereas the superburst itself does not dif- 
fer substantially from models with a carbon-rich atmo- 
sphere. After the superburst all burning in the atmo- 
sphere is stable, and bursts are quenched until the en- 
velope has cooled down from the superburst. At that 
point burning first becomes marginally stable, followed 
by weak bursts with short recurrence times. Over time, 
recurrence times lengthen, and the burst properties re- 
turn to those from before the superburst. 

4.1. Superburst 

IKeek fc HegeH ()2011[ ) found for a model with the 
same accretion rate and a slightly higher Qb — 
0.13MeVnucleon~^ a recurrence time of 1.70 years, 33 % 
longer than the recurrence time found in the present 
study. With Qb = 0.1 MeVnucleon"^ one would expect 
our cooler model to produce a somewhat longer recur- 
rence time. For the models in the current paper, we use 
a different implementation of accretion and, most impor- 
tantly, of the compressional heating from the accreted 
material. A further difference is the prescription used 
for electron conductivity. 

The different phases di splayed by the superb urst are 
the same as described by IKeek fc Hegen ()201in . includ- 
ing the shock breakout, the precursor, and the two- 
component power-law decay. A difference lies in the 
height of the shock breakout peak and the oscillations 
at the start of the precursor, where super-Eddington 
luminosities are reached. The peak luminosity has a 
strong dependence on the resolution at the surface of 
the model (Keek & Heger 2011). The bottom of the 
outer zone of our models is at a column depth of y ~ 
lO^gcm^^, whereas the photosphere of a neutron star 
is typically located at y ~ Igcm^^, so it is likely that 
our simulations have not resolved this fully. A previ- 
ous study also finds the shock breakout to be super- 
Eddington, but likewise docs not fully resolve the photo- 
sphere (Weinber g fc Bildsten 20071) . This underlines the 
importance of the treatment of the outer layers during 
these hydrodynamic events. Whereas our non-relativistic 
simulations employ diffusive radiation transport, a rela- 
tivistic hydrodynamic model that includes full radiation 
transport will be much better suited to simulate these 
processes in detail. Such a model will also be able to 
determine whether any material is lost from the neutron 
star during the super-Eddington luminosity phases. 

4.2. Precursor 



After the initial oscillations, the precursor luminosity 
settles at the Eddington limit, £Edd, producing a plateau 
in the light curve. The height of the plateau is different 
for the three atmosphere compositions because of the de- 
pendence of iEdd on the opacity, k: Lsdd oc k^^. The 
opacity in the neutron star photosphere is typically as- 
sumed to be dominated by Thomson scattering, which 
depends exclusively on the hydrogen mass-fraction. Both 
models He and C are devoid of hydrogen, but they reach 
the Eddington limit at slightly different luminosities. 
The opacity in the outer zones of the models continue 
to increase toward the surface. It is likely that the true 
photospheric values are not yet reached, which may ex- 
plain the variations in ^Edd- 

The precursor in model C is predominantly powered 
by the fallback of shock-heated material. The fluence 
of this precursor is about half that of the precursors in 
the models He and H, indicating that even in those cases 
fallback has a substantial contribution to the precursor, 
with hydrogen and helium burning accounting for the 
rest of the fluence. 

The superburst ignited at a relatively early phase in 
the cycle of hydrogen/helium bursts. If it had ignited 
at a later phase, the respective helium-rich and so- 
lar composition columns would have been several times 
larger, allowing for somewhat longer precursor bursts. 
Note that the brightness of the precursors would not 
change, as it is set by the Eddington limit for the 
atmosphere composition. Precursor bursts that reach 
the Eddington limit with durations of several seconds 
have been observed in the few instances when the start 
of the superburst was observed and the data were of 
sufficient qualitv fe.g.. IStrohmaver fc Markwardtl [20021 ; 
iStrohmaver fc Brownll2002f ). 

Previous models of a superburst below a helium atmo- 
sphere skipped the fallback , losing an important energ y 
source for the precursor ([Weinberg fc BildstenI l2007f ). 
The shock breakout left a flat temperature profile, de- 
laying the burst from helium burning by a thermal 
time scale of several seconds. The precursor of model 
He arrives on a dynamical time scale of ~ 10~^ s, 
and is powered by both the fallback and the subse- 
quent thermonuclear burn ing of helium. Furthermore, 
IWeinberg fc Bildstenll2007l argue that the thermonuclear 
runaway of helium is triggered either by the shock if the 
helium layer is sufficiently thick, or by the carbon de- 
flagration reaching the helium layer. In our simulations 
the carbon deflagration does not reach that far out, and 
the shock does not trigger the helium burst immediately. 
It is rather the fallback that heats the atmosphere suf- 
flciently to burn helium, even if the helium layer is thin 
(Fig.®. 

4.3. Burst Quenching 

The superburst heats the atmosphere such that all sub- 
sequently accreted hydrogen and helium burn stably, and 
no bursts are produced. Burst quenching continues until 
the envelope cools down sufficiently for burning to be- 
come unstable again. Using the analytic estimate from 
iCumming fc MacbethI (|2004[ ). one predicts for the condi- 
tions of our model a quenching time of 5.5 days for model 
H and 2.2 days for model He. Here it is assumed that 
3a is solely responsible for the stability of thermonu- 
clear burning. Model H, however, yields the shortest 
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quenching time of 1.1 days, whereas in model He bursts 
are quenched for 11.3 days. The latter is over five times 
longer than the analytic estimate, which may indicate 
that the approximations employed in the estimate are too 
strong. In contrast, we find a five times shorter quench- 
ing time for model H. In this case 3a is not the sole 
process responsible for the resumption of bursts, but the 
CNO breakout reaction ^^0(a,7)^^Ne plays a key role. 

The importance of the ^^0(a, 7)^^Ne reaction for regu- 
lating the CNO abundance i n the co ntext of X-ray bursts 
was stressed before by F isker et al.l (|2006.) , and indicated 
as an im portant factor in the sta bility of the burning 
processes (jFisker et al.ll2007l [20081 ). Our study shows its 
importance in determining the duration of the quenching 
period. 

Mixing due to rotation, or a rotationally induced mag- 
netic field, can also influence the stability of thermonu - 
clear burning ()Piro fc BildstenI [2007t [Keek et all [20091 ). 
Rotation was not taken into account in this study, but 
could lead to a somewhat longer quenching period. 

For the ignition column depth of our simulations, y = 
1.1 X lO^^gcm-^. lKeek fc Hegeii (|2011f ) showed that the 
fraction of the superburst energy that is lost in neutrinos 
is small, whereas for ignition at larger depths neutrino 
losses are substantial, and the photon fluence reaches 
a maximum value. This maximum fluence implies that 
there is a maximum burst quenching time, which is longer 
than the values found in our simulations. 

4.4. Marginally Stable Burning 

When burst quenching ends, burning becomes 
marginally stable, producing oscillations in the light 
curve. This has been observed from a small num- 
ber of sources as mHz quasi-p e riodic oscillations (mHz 
OPOs: iRevniytsev et al] 120011 ; lAltamirano et al.l I2008L 
iLinares et al.l l2012f ). and is associated with a burning 
mode at the transit ion of stable and unstable burning 
(jHeger et al.ll2007H ). Using the one-zone analytic ap- 
proximation from Equation (11) in lHeger et al.l (|2007bD . 
and substituting the values appropriate for model He 
= 594keVnucleon-i, Tg = 0.54, yg ^ 1-0, and 
m = 0.3 mEdd), one obtains the oscillation period Pose — 
25 minutes, which is close to the 20 minute period in our 
simulations. Marginally stable burning in model H has 
a period of 5 minutes, which is close to Pose — 6 minutes 
predicted by the analytic approximation with Tg — 4.2, 
2/g = 0.60, and r h = 0.3 mEdri- The analytic estimate 
from iHeger et al.l ()2007bl ) considers only accretion, ra- 
diative cooling, and a single thermonuclear process, the 
rate of which has a similar temperature dependence as 
the cooling rate. In the case of model He, where 3a is 
the single dominant thermonuclear process, this approxi- 
mation seems valid. It is interesting that it also provides 
a good estimate of Pose in model H, where the nuclear 
reactions proceed through a complex interplay between 
the 3a, hot-CNO, and CNO breakout processes. 

4.5. Burst Resumption 

As the envelope continues to cool, the oscillatory 
burning transitions into bursts. At first their flu- 
ence is relatively low and the recurrences times are 
short. Over time (almost four months for model He; 
one month for H) the burst properties regain their 



pre-superburst values. The light curves of the later 
bursts, as well as of the pre-superburst bursts, are qual- 
itatively similar to tho se reported from recent multi- 
zone simulations (e.g., 
[200l [Fisker et aLl l""^ 



Joss &: Lil [198(1 iWooslev et all 
Jose et ani2010l ). w hich provide 



very good agreement with observed bursts (IHeger et al.l 
l2007al:[h7VZand et al.l[200l ICvburt etamom . 

The transition from stable burning, to bursts, via 
marginally stable burning (mHz QPOs), and weak 
bursts, has not been observed directly after a su- 
perburst. In this respect the bursting behavior of 
ICR J 17480 -2446 in the globu la r cluster Terzan 5 is 
interesting (iMotta et all |201H lLina.res et all 1201 Ibl; 
Chakrabortv fc Bhattacharvval 120111 ; ILinares et al.l 



20121) . During an outburst in late 2010, its accretion 



rate varied by more than a factor five, and it displayed a 
continuous transition from bright bursts to weak bursts, 
to mHz QPOs, and back. The changes in the burst 
properties are quite similar to those in model H. The 
weak bursts have a longer rise and recurrence times as 
short as a few minutes. The bursts become brighter with 
longer recurrence times and longer decay profiles. The 
decrease of the accretion rate causes the atmosphere to 
cool, similar to the cooling after the superburst in the 
models, but it also changes the burst fuel accumulation 
time. This is more complex than the situation in the 
simulations, and we cannot compare the time scale on 
which the changes in burst behavior take place, as it is 
set by the variations in the mass accretion rate, instead 
of the cooling time scale of the envelope. 

4.6. Carbon Production and Destruction 

One of the biggest challenges in superburst theory is 
the creation of the correct amount of carbon. Fits to 
superburst light curves det ermine the carbon mass frac- 
tion to be close to 20% (jCumming fc Macbeth! 120041; 
iCumming et al.l [20061) . which is the amount we adopted 
for our models. Models of helium bursts, or mixed hy- 
drogen /heliu m bursts, however, produce less than half 
of that (e.g.. IWooslev e t al. 2004). It is speculated that 
both bursts and stable burning of hydrogen and helium 
are required to produce enough carbon for superbursts. 
In fact, most superbursting sources exhibit this c ombi- 
nation of burning behavior (jin 't Zand et al.ll2003D . 

Our models exhibit both burning regimes, and the next 
superburst will ignite at the bottom of the layer that is 
accreted just after the simulated superburst. Because of 
the relatively high temperature in the envelope after the 
superburst, however, most of the produced carbon burns 
once it is compressed to higher densities over time. Sim- 
ulations over a longer time are required to determine how 
much carbon can survive when the envelope cools down 
further, and the carbon fusion reaction rate is reduced. 
Model He might produce a carbon-rich layer where the 
next superburst can ignite, whereas in model H most 
c arbon is destroyed. 

iMedin fc CummingI ()2011[) recently suggested that a 
larger carbon fraction can be obtained at the superburst 
ignition depth by separating out the carbon from heavier 
elements by freezing at the crust-ocean interface. 

4.7. Observations after Superbursts 

The detection of a burst is a certain upper limit to the 
quenching time. Table [2] presents the current observa- 
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tional upper hmits on burst quenching for ah known su- 
perbursts, based on the first detection of a burst fohowing 
each superburst. We identify an upper Umit of 15.0 days 
for 4U 1636-536 using the MINBAR burst catalog. This 
is the shortest reported value apart from the 2.2 day limit 
for GX 17+2. which accretes at a much higher rate than 
all other superbursters. 4U 1636-536 accretes hydrogen- 
rich material at a rate of approximately 0.1 MEdd- The 
15.0 day upper limit is well above the 1.1 day quench- 
ing time we find for model H, which has a three time 
higher mass accretion rate. The shortest upper limit 
for UCXBs, which may ac crete helium-rich ma terial, is 
18.6 days for 4U 0614+091 (|Kuulkers et alllMoh . Model 
He has a lower quenching time of tquench = 11.3 days, but 
has a 30 times larger mass accretion rate. 
Gumming & Macbeth (2 004) estimate the dependence 

of iqucnch ou the accretiou rate to be tqucnch oc M^^/^. 
Using this to correct for the differences in accretion rate 
between the models and the observations, model H's 
^quench — 6.6days is still below the observational up- 
per limit, but model He's iqucnch — 28.9 days lies above 
the upper limit. Gonsidering, however, that the analytic 
approximation does not reproduce the results from our 
multi-zone models well, the accretion rate dependence 
here employed requires further scrutiny before conclu- 
sions can be drawn. 

The two bursts following the superburst from 
4U 0614+091 have the same peak flux and suggest a 
trend towards longer bursts at later times, just as pre- 
dicted by model He. 

For the hydrogen-accreting sources we combine the ob- 
servations following superbursts, that are reported in 
MINBAR. There is a strong indication for burst quench- 
ing during the first 4 days after superbursts. A burst 
observed after 26 days firmly constrains the tj-ccur at a 
pre-superburst value. Furthermore, the lack in variation 
of the properties of bursts observed from 4U 1636-536 
suggests that pre-superburst bursting behavior had re- 
turned already after 15 days, which is shorter thaii the 
35 days predicted by model H. lGumming et all (|2006D de- 
rive from fits to the light curve of a different superburst 
from the same source an ignition column depth that is 
lower by a factor 0.45 compared to the ignition depth in 
our models. This implies a shorter cooling time scale af- 
ter the superburst (tcooi oc Ihln^ iGumming fc MacbethI 
l2004f ). and hence a faster return to normal bursts of 
19.3 days, which is closer to what we infer from the ob- 
servations. 

For all but a few days, we derive frecur > 1 hr. This 
leaves little room for the short recurrence times pre- 
dicted by our models and observed from IGR J17480- 
2446. The majority of the observations during the first 
week after superbursts were, however, performed with 
the BeppoSAX WFGs, which did not have enough sensi- 
tivit y to detect o scillatory behavior or weak bursts (see, 
e.g.. iKeek et"an[20Tol for a comparison of the detection 
of weak bursts with the RXTE PGA and BeppoSAX 
WFGs). Alternatively, the bursting behavior after a su- 
perburst may resemble that close to the transition to sta- 
ble burning ob served at higher M, w hen the burst rate 
decreases (e.g.. lGornelisse et al.ll2003l ). 

Both models H and He predict a period when weak and 
brighter bursts alternate. This behavior has not been 



observed, which is not surprising considering the short 
time that it is expected to take place. For model H this 
burning mode lasts a few hours, whereas for model He it 
continues for 2.3 days. It is, therefore, easiest to observe 
this behavior from a UGXB. 

Because observations after a superburst are often few 
and far apart (Fig. we have to combine the data 
for superbursts which differ by up to a factor 6 in j/jgn 
(IGumming et al.l [20060 . and have corresponding diifer- 
ent quenching and cooling times. A campaign of frequent 
and longer observations during the month following a su- 
perburst would, therefore, be very important in improv- 
ing our current inability to accurately predict changes 
in bursting behavior. Now that RXTE has ceased its 
operations, there is no observatory available with a suf- 
ficiently large effective area and time resolution to study 
the mHz QPOs and weak bursts, as well as the brief su- 
perburst precursors. A future mission such as LOFT, 
that greatly improves on collecting area, will be able to 
make important steps forward in our unde rstanding of 
therm onuclear processes on neutron stars (|Feroci et al.l 
l20ll . 

5. CONCLUSIONS 

We create one-dimensional multi-zone models of a 
neutron star envelope undergoing a superburst (carbon 
flash), in the presence of an atmosphere of either pure 
helium or of solar composition. The latter is the first 
model of its kind. After the superburst we continue the 
simulations to study burst quenching as well as the re- 
turn of normal Type I X-ray bursts, simulating over 900 
hydrogen or helium hashes per model. 

The heating of the atmosphere by fallback of shocked 
material generates a precursor burst that reaches the Ed- 
dington limit. Any available hydrogen or helium ignites, 
which extends the duration of the precursor. 

After the superburst, the atmosphere is sufficiently hot 
for bursts to be quenched, and all hydrogen and helium 
to burn in a stable manner. In a pure helium layer mostly 
carbon is produced, and in a layer of solar composition 
stable rp-process burning creates mostly germanium and 
gallium. As the envelope cools down from the superburst, 
bursts reappear in the light curve. In the helium atmo- 
sphere this happens after 11.3 days, when the 3Q;-process 
becomes unstable. We find that in a solar composition 
quenching ends much sooner, after only 1.1 days. In this 
case it is the interplay between 3a, the hot GNO cycle, 
and the ^^0(a,7)^^Ne breakout reaction that leads to 
unstable burning. 

In both models, at the transition from stable burn- 
ing to unstable burning, oscillations are produced in the 
light c urve due to marginally stable burning (Hcgcr et al.l 
2007b) , followed for a brief time by alternating weak and 
brighter bursts, which have longer and shorter recurrence 
times, respectively. The latter bursting mode has not 
been observed yet. Eventually the weak bursts disap- 
pear. In the helium atmosphere the bright bursts are 
immediately as bright as those before the superburst, 
whereas in the solar atmosphere burst peak luminosities 
grow as the envelope cools further. For a few days left- 
over helium above the burst ignition column depth burns, 
pausing the cooling. Afterward the burst durations and 
recurrence times increase back to the pre-superburst val- 
ues over the course of one month for solar composition, 
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and four months for helium composition. 

The transition from burst quenching to bursts after 
a superburst has not been directly observed. Using the 
MINBAR catalog we identify the shortest reported upper 
limit to the quenching time of 15.0 days for 4U 1636- 
536 (with the exception of GX 17+2), and we derive 
further constraints on the time scales for quenching and 
the return of bursts. The short recurrence times found 
by the simulations are disfavored, but not excluded. The 
transition between the different burning regimes that we 
describe exhibits, however, strong similarities with bursts 
observed from the transient burste r IGR J17480-2446 in 
Terzan 5 (e.g.. iLinares 

The authors thank K. Chen for helpful discussions. 
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